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Abstract 

We investigate the pion photoproduction off the nucleon in the A region in the 
framework of an effective chiral Lagrangian including pions, nucleon, and A(1232). 
We work to third order in a small scale expansion with both m n and Ma — Mw treated 
as light scales. We note that in the A region, straightforward power counting breaks 
down as the amplitude becomes very large. To deal with this problem, we suggest that 
the appropriate way to compare the theoretical calculations with experimental data is 
via weighted integrals of the amplitudes through the A region. 
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1 Introduction 



The A(1232) resonance enjoys a special status in the family of the nucleon resonances. It 
lies only about 300 MeV above the nucleon ground state, nicely isolated from the plethora of 
other densely populated resonances at higher energies, and only strongly couples to the Ntt 
system. Therefore, the matrix of physical current-like (N\J^ m \A) is supposed to be easily 
extracted from experimental observables and expected to shed light on their structures. 

For example, in simple constituent quark models, the interaction between quarks is be- 
lieved to have a tensor part due to the color hyperfme interaction [|l|. One consequence is 
the admixture of higher orbital angular momentum components in s-state quark wave func- 
tions of nucleon ground state and A. The d-state mixture allows for an electric quadrupole 
E2 transition in 'yN — > A(1232) excitation, which is otherwise a pure magnetic dipole Ml 
transition. The Ml and E2 transitions can be directly excited by photons and the subse- 
quent pion decay can be observed. The amplitudes in 7rN final states are usually denoted 
by E( and M( where E and M are respectively electric and magnetic multipoles, I is the 
orbital angular momentum of the pion, the ± signs refer to the total angular momentum 
j = I ± 1/2, and I is the isospin of the 7rN system. Therefore, in this class of models, the 
ratio Rem = E2/MI = E^/Mf^ is related to the tensor component in the interaction of 
quarks. The simplest noninteracting SU(6) quark model predicts Rem=®', the nonrelativis- 
tic constituent quark model with gluon exchange predicts a rather small and negative result 
ranging between -0.008% to -2%[|],@. 

Many baryon models have the pion cloud playing an important role. Because of the 
derivative coupling pions are relatively efficient at generating strength for the E2. Cloudy 
bag models lead to —2% < Rem < 00. Larger negative values are given by Skyrme 
and other hedgehog models:— 5.9% < Rem < — 2.5%|], and exchange current yields: 
Rem = — 3.5%0. Thus Rem becomes a sensitive test for different models for baryon 
structure. 

On the other hand, for several different reasons it is important to understand the N — > A 
EM transition and more generally, the nucleon and A-isobars, without reference to any 
particular baryon model. For example, to study the photoproduction of the mesons off 
complex nuclei ||, it is necessary to use single nucleon information such as '-/AN vertex 
strength as input. Such investigations help to clarify how the effective degrees of freedom, 
nucleons, mesons, and A-isobars play their roles in nuclei. A-isobars by themselves are 
particularly interesting objects in the large N c QCD because they will be degenerate with 
the nucleon if N c exactly goes to infinity. By a simple argument based on N c counting rules 
applied to the tcN scattering, it was shown that large N c QCD with only I=J=l/2 nucleon 



interacting with pions is inconsistent [10]; other states in the tower of I=J=l/2,3/2 ... must 
be included to satisfy the consistency conditions. Therefore, the study of N — ► A transitions 
may help provide insights into the large N c QCD. Moreover, Rem was predicted to be unity 
in the domain of perturbative QCD||; thus the study of N —>■ A EM transition may provide 
a window into the breakdown of pQCD as momentum transfer drops. 

3 /2 

It is difficult to determine the resonant E 1+ for two reasons. From the experimental 
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side, the precise measurement of this quadrupole amplitude is difficult due to its smallness 
compared to dominant magnetic dipole amplitude. Recently there has been substantial 
progress so that precise measurements of spin observables are now possible. At Mainz [11] a 
p(^p)7r° measurement was performed with tagged linely polarized photons produced at 855 
MeV Mainz Microtron MAMI; the differential cross section and photon asymmetry £ were 
measured simultaneously for the pit and mr + channels. They took Rem = ImE\+ / ImM^ 
at the peak of the A(1232) resonance, and the reported value is Rem = —(2.5 ±0.2 ±0.2)%. 
There are also analyses of other groups based on their and BNL's data[l2"|-[I6|. 



However, there is a more serious problem on the theoretical side: The non-resonant 
contribution in E±+ cannot be separated directly by the measurement fi"8|| . Thus, the direct 



comparison between the experimental Rem with the results of the calculations based on the 
specific baryon models are meaningless unless one extends the model to include continuum 
states in a consistent way. Various methods have been proposed for the isolation of the 
resonant part from the measured multipoles results but their answers vary greatly. Different 
prescriptions result in different kinds of definitions of the resonant E2/M1 ratio [jHJ . 

For example, in the effective Lagrangian approach one main difficulty of this separation 
is due to the unitarization. The Born terms and the A excitation treated as the leading tree 



graphs, do not fulfill the requirement of unitarity[20]. The unitarity may be put in by hand 



but different unitarization prescriptions lead to different separations between the resonant 



part and background contributions [EH]— [B5| 



On the other hand, models ||24||,||25||. which treat pion photoproduction dynamically (i.e. 



solving the corresponding Lippman-Schwinger equation for a given irN interaction) are au- 
tomatically unitary since the rescattering process is included. Such models also provide a 
basis for the analysis of the role of final state interactions. However, to solve the dynamical 
equation, some phenomenological form factors must be included. They are needed to regu- 
larize the driving terms of the interaction. Since these form factors are put in by hand, they 
also make the separation between the resonant part and background contributions model 
dependent [f26| . 



Recently, the "speed plot technique" was also applied to this problem [[17]]. The "speed" 
SP of the scattering amplitude T is defined by: 

Here W is the total cm energy. In the vicinity of the resonance pole the energy dependence 
of the full amplitude T = Tr ± Tb is determined by the resonance contribution: 

T R {W) 



M R -W-iT R /2 



while Tb is a smooth function of energy. Application of this method to the amplitude of 
Tiator et al. derived by fixing the t dispersion relation gives 



R A = . , = -0.035 - 0.046z. 
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Th. Wilbois et al.[0 suggested a more sophisticated way to implement this and found Ra = 
—0.040 — 0.0472. Note this is a complex and energy-independent quantity; unfortunately 
it remains unclear how Ra could be compared to any microscopic baryon model, and it is 
difficult to determine what, if anything, about the chiral dynamics of ttAN system. 

Some of the ambiguities mentioned above will be ameliorated by using an approach based 
on chiral perturbation theory. Unitarity is guaranteed, at least perturbatively, because loop 
graphs are included. Since dimensional regularization can be implemented straightforwardly 
in this approach, one avoids the model dependence inherent in the introduction of phe- 
nomenological form factors. Therefore ChPT provides a model-independent picture of pion 
photoproduction — at least to the extent the expansion converges. Many previous ChPT 
calculations were limited to the threshold region because A-isobars were not included ex- 
plicitly. In such calculations, the A effects are supposed to be included in the form of local 
counterterms. The essential physical idea is that the delta propagator is treated as though it 
was shrunken to a point and the energy dependence of A-isobars are reproduced by higher 
dimension operators which are suppressed by 1/(Ma — Mat) [PTfl , P8| . Of course, such a 



framework cannot be used for our problem since we wish to work in the A region. 

It has been recognized, even in the threshold region, that the A is a low energy ex- 
citation and it is presumably sensible to include it dynamically, creating a more general 
effective field theory than simple ChPT. This is the spirit of the original work of Jenkins and 
Manohar p9| , [pOj] , and the applications by Butler, Savage and Springer to the SU(3)xSU(3) 



case 



31], [3^| although their works are never beyond the leading order. 



Recently, Hermmert et al.|p3|-||35| have formalized such an approach. They have de- 
veloped a consistent power expansion scheme, the so-called "small scale expansion", which 
allows for nucleon and A-isobar degrees of freedom to be treated simultaneously in an SU(2) 
effective chiral Lagrangian. Whereas in conventional HBChPT one expands in power of 
external momenta in analogy to the meson sector; a phenomenological expansion was set 
up in the small scale e. "Small scale" denotes the soft momentum, the pion mass or the 
mass difference A = Ma — M^. One natural reason to do this is: A now is treated as a 
new dimensionful parameter which stays finite in chiral limit but is nevertheless of compa- 
rable size to m n in the real world. To assert the accuracy of this novel approach one has to 
systematically calculate observables and compare the resulting predictions. The 7iV — > A 
transition seems to be a promising case to implement such a scheme. [f49| 

There is a potential problem, however. The perturbative power counting scheme fails 
in producing calculations for momenta transfers in the A region. This is unfortunate since 
this is precisely the region where we wish to work. The reason for this failure is quite 
obvious. The generic power counting has the A propagator behaving as 0(l/e), where e 
is the small dimensionless parameter, e = m^/A, A/A. For generic low momenta this is 
valid. As one approaches the A resonance, the propagator, treated at lowest order diverges, 
spoiling the power counting. One may hope to cure this by including the A self-energy in 
the propagator. The imaginary part of the self energy remains finite as there is a physical 
decay channel. There is a difficulty with this approach; namely, certain graphs are iterated 
to all orders — the A self energy insertions — while others are not. This makes it very difficult 
to assure that one has a systematic power counting scheme. However, even if the propagator 
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is anomalously large in the neighborhood of the resonance, we note that integrals of the 
propagator times smooth functions over regions of order eA are order 1 / e — the same as the 
generic A propagator. Clearly, this means that so long as we are only sensitive to integrals 
over the propagator the power counting is safe. Loop diagrams with A's as intermediate 
states are of this type. However, there are also contributions coming from graphs such as 
the Born graphs, in which the four momentum flowing through the propagator is fixed by 
the external momenta. When those external momenta are such that four momentum in the 
A propagator are close to the pole the power counting has broken down. 

We propose a possible way to avoid this problem. The power counting scheme used simply 
does not allow us to accurately calculate the amplitudes in the vicinity of the A pole. How- 
ever, if we limit our ambitions to asking questions about integrals of the amplitudes through 
the A region the power counting scheme remains viable. Accordingly, our proposal is that 
one should not directly compare calculated amplitudes with experiment. Rather one should 
extract the amplitudes from experiment, estimate weighted integrals over the amplitudes 
and compare these integrals of the experimental amplitudes with theoretical ones. In this 
way, we can make predictions of quantities in the A region — albeit integrated quantities — for 
which the power counting scheme is viable. We note that this approach has some important 
limitations. The most obvious one, of course is that we cannot make a direct prediction of 
the experimental observables. There is also an important practical limitation. We do not 
make predictions for integrals of differential cross sections but for integrals of amplitudes. 
In order to do this, one must fix the amplitudes from the measurements. Unfortunately, the 
various spin-dependent differential cross sections each depend on several amplitudes. One 
needs to make several independent measurements to extract the amplitudes |p4||. To the best 



of our knowledge none of the amplitudes has been extracted from the experimental quantities 
to date. This means we can not presently use the methods discussed here to compare with 
experiment. However, future spin-dependent measurements could alter this situation. 

This paper is organized as following: In Sec. 2 the formalism of HBChPT and the steps 
to include the A degree of freedom are briefly sketched. Sections 3 and 4 describe the loop 
and Born graphs. Renormalization is discussed in section 5. Finally, in Sec. 6, the method 
for comparing the theoretical calculation with experimental data is reported; some related 
issues and further prospects are also discussed. 



2 Effective Lagrangian 



The starting point for our approach is the most general chiral invariant Lagrangian involving 
relativistic spin 1/2 and spin 3/2 fields: 

C = C N + C A + (C AN + h.c). (1) 

Although the expressions of the standard conventions of SU(2) HBChPT exists in literature 
as in Ref. [34], we believe that it is useful in establishing our notation to give relevant 
expressions here: 
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(17) 



When the only external fields are photons, simplifies to: 

f% = e(d„A v - d v AJ{uQu^ ± vtQu), (18) 

Q=diag(l,0). denotes external vector field and a M denotes axial field. Here x{ x ) — 
s(x) +ip(x) includes the explicit chiral symmetry breaking through the small current quark 
masses, s(x) = B diag(m u ,mA B = with F n the weak pion decay constant. Here 

isospin invariance (m u = m d ) is assumed. 

The spin-3/2 field is represented via a Rarita-Schwinger spinor[3(|, i.e, a vector-spinor 
field ty(j,(x) satisfies the equation of motion: 

(i lv d v - M A )^(x) = 0, (19) 

with the subsidiary conditions: 

!^(x) = 0. (20) 

This subsidiary condition ensures that the physical particle created by the field is spin-3/2 
as opposed to 1/2. Following Hermmert et al.[]33|| the Lagrangian of A-isobars is given by: 

C A = #A«^, (21) 
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+ f ^ ^75 + f (7m ^7,)75 + f 7m ^757,]- 



(22) 



£at is well known in HBChPT, to the third order it is given by H38| , [j39|l . Only nine of the 
terms in are relevant in our problem. As Fettes et al.jS^] pointed out that there is one 
counterterm too many in in p8[, however the extra term 0± is irrelevant in the pion 



photoproduction process, so we still adopt the basis of [p8 |. 



The interaction of the A field with photons, pions, and nucleons is given by the following 
effective Lagrangian [ H | , [ ¥2" ] : 



AN 



r (l) , n 



(2) 



c 



(3) 



c 



(2) 

txAN 



m N 



(23) 
(24) 
(25) 



We need all AiV vertices is up to 0(e 2 ). However the following term must be included 
even though it is apparently third order, because it produces one C(e 2 ) 7rAiV vertex in heavy 
baryon limit: 



C 



(3) 
:AN 



ivrAiV T 



m 



N 



(26) 



The 7 AiV vertices are given by following terms: 



A 1 ) 

l ~~ i AN 
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2m 
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4>i®^( x oh P j5f+ u i>N + h.i 
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(2m N ) 



-M < d t iu{y2)^\ P lbf+ P ' U i , N + h.c 



^0^(^3)^75/+ U '' P iPn + h.c 



2m 2 N 



(27) 
(28) 
(29) 

(30) 

(31) 

(32) 



Again, we have to include the following terms which were not be considered in references 
such as fRI , because they produce 0(e 2 ) 7 AiV vertices: 
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'"jAN ~ ""jAN 7 A AT) \ 00 J 



4 3 A a i = ^W®»M lxl5 fi» x > p D p Tp N + h.c. ( 3 I ) 



£§Sr = 7^^i^^{ti)l\l^f+ P,VD P ^N + h.c. (35) 



The tensor Q^ u (z ) = g^ u + z ol^lu was introduced by Peccei|43[] as the most general form 
obeys the invariance under the contact transformation. The so-called "off-shell" parameters, 
z and xq, yo, etc. govern the couplings of the off-shell spin- 1/2 fields to external fields. There 



has been considerable controversy in finding conditions to fix these parameters ||43|| - |4"5j| . We 



find however, that at the end of our calculation that final result is completely insensitive to 
the off-shell parameters. This is quite reasonable. After all, as a general rule in a consistent 
scheme, off-shell effects are not observables; while they enter calculations in intermediate 
stages they should be canceled in the final answers. This is clear in the case of A from the 
KOS theorem [Tf \ and presumably applies for the other off-shell parameters as well. In this 



context it is useful to note that Tang and Ellis [46] explicitly showed such independence for 



a somewhat simpler system than the one considered here. 

The next step in this approach consists of identifying the "light" and "heavy" degrees of 
freedom of spin-3/2 fields. The procedure is in analogy to the case of spin- 1/2 fields. The 
situation becomes a little more complicated here because of the off-shell spin-1/2 components 
associated with the Rarita-Schwinger field. One can identify the "light" component as: 

T; = P v + P^ v W(x) exp(iM N v ■ x). (36) 

The remaining component: 

Gl = (g^ - P+P*£ )lx M(x) eM*M N v ■ x) (37) 

The nucleon field is also split as "light" and "heavy" components as in the usual HBChPT: 

N(x) = P+ip N exp(iM N v ■ x), (38) 
h(x) = P^ipN exp(zMivt> • x). (39) 

The general Lagrangian now take the form: 

L N = NA N N + (hB N N + h.c.) - hC N h, (40) 

L AN = TA AN N + GB AN N + hD NA T + hC NA G + h.c, (41) 
L A =fA A T+(GB A T + h.c.)-GC A G. (42) 
The matrix A^, B A n..., admits a small energy scale expansion of the form: 

A N = A$+A% ) +AP + , (43) 



S 



where A$ is of order e n . As emphasized in the introduction, we denote by e small quantities 
of order p, like m n or the soft momenta, as well as mass difference A = M A — M n . This 
mass difference is distinct from the pion mass in the sense that it stays finite in the chiral 
limit but vanishes in the large N c limit. However, in the physical world, A and m n are of the 
same scale — differing by only a factor of ~ 2. We therefore use a simultaneous expansion 
in both quantities. It is only through this small scale expansion that we may obtain a low 
energy expansion of 7rAiV system. Such an expansion is in the spirit of large N c ChPT since 

The standard procedure is to integrate the heavy components. This results in a relatively 
simple effective action: 

s eff = J d 4 x(TA A T + NA N N + [TA AN N + h.c.)), (44) 



with 



Here 



A A = A A + -f D NA -f C N l D NA + 7o5 A 7oCX 1 5 A , (45) 

A N = A N + l0 B ] Nl0 C N l B N + 7o£ A 7v7oC A l5 Aiv, (46) 

A AN = A AN + l0 D j NAl0 C^ l B N + 7o£ A 7oC A 1 J B AJV . (47) 

C N = C N — Ci VA C A 1 7oC , j VA 7o, (48) 

B N = B N + C NA C A X B AN , (49) 

Dna = D NA + C NA C^B A . (50) 



Note C A is a 5 x 5 matrix pi. The effect of heavy degrees of freedom shows up at order e . 



In order to calculate a given process to order e n , it is sufficient to construct matrix A to the 
same order, e n , B and D to order e" -1 , and C to order e n_2 . 



3 Loop Graphs 

To order e 3 , only one- loop graphs with the lowest order of vertex need be considered. The 



vertices we need are from A K N ,A A and A AN . 



> m ' : a,„> » |J| 



Atf=iv-D + g A S ■ u, (51) 

A { AN = g V ANV>l, (52) 

A$ = -[iv D ij - AS ij + g x S ■ u l3 \ 9iiv . (53) 

Here = ^jsv^v" denotes the Pauli-Lubanski spin vector. From A$, A A we determine 
the propagators in momentum space with residual soft momentum: r li = p^ — M^v^. 

S 1 / 2 ( v . r ) = L_ (54) 
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(55) 



with denoting the spin 3/2 projector in d- dimensions: 



4 



(56) 



d-4 



and 




(57) 



denotes an isospin 3/2 projector. At present, the coupling constant g n AN has not been 
extracted from data within the context of a systematic implementation of the small scale 
expansion. 

For simplicity, we work in the center of mass. We choose v M = (1, 0, 0, 0) which is 
equivalent to working in (or close to) the center of mass frame, and the entire calculation 
is done in Coulomb Gauge: e • v—0. This choice greatly reduces the number of graphs 
which contribute. The disadvantage is that the gauge invariance is no longer manifest. The 
properties of the "light" components of the delta: 



also ensures that many graphs vanish identically. Note that imAN vertex, unlike ttttNN 
vertex, starts from 0(e 3 ), and since A-isobars must be the intermediate states, the 7T7rAA 
vertex and -kit AN vertex do not enter our calculation. In Figs. 1, 2 all graphs are shown 
(l-a,b,c are grouped in gauge invariant classes). 



4 Born Graphs 

The Born graphs contributing to third order are shown in Figs. 3, 4 with the vertices given 
in Table 1: 



ltl T^{x) = v^\x) = 0, 



and 



SpT^ix) = 0, 
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Table 1: Vertices for Born graphs in pion photoproduction. 



The 7 iVA vertex is the focus of our study here. The leading order 7 iVA vertex is a 
pure Ml one, but through 1/M expansion, the coupling constants G\ and G2 still contribute 
to E2 transitions. The 0(e 2 ) vertex with G4 is the same with the one with G5 and G7, 
therefore we use G4 to represent the G4 + G5 + G-j. Similar situation occurs in the 7rAiV 
vertices: the C(e 2 ) vertex with K^an is the same with g^AN, so we also use g n/ \N to represent 
the tin an + 9-kAn- Of course if we go to higher order calculation, these different coefficients 
should contribute in different way. 

Note also that the heavy and off-shell spin 1/2 component of A-isobars modify the NN 
vertices, up to third order giving nonvanishing contributions. The off-shell parameters show 
up through the 1/M expansion. The combinations hi = y§G<i — 2xqG\ — 2Gi, and h 2 = 
—2xqGi + 110G2 appear in the 0(e 2 )7AiV vertex. However, the amplitudes depend only on 
hi — h%, which is 2Gi and is independent of xq, yo. (Actually the vertex itself can be shown 
as independent of xq, yo if the on-shell constraint: S^T- 1 = is implemented, then our vertex 



is identical with the one in [49] 



The off-shell parameters are also encountered in the 0{e 2, ) r )TrNN vertex, but they cannot 
be distinguished from 6 9 ; thus, in our calculation, the values of off-shell parameters are irrel- 
evant. As discussed earlier, this is expected on very general grounds. Also the counterterm 
in £(*) appears in Fig. 3-C-7. 
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5 Renormalization 



The one-loop diagrams contain divergences and have to be regularized and renormalized. 
Here we will use dimension regularization. The unrenormalized coefficients are then related 
to the renormalized (scale dependent) ones according to 

h = b r M + (4vr) 2 AL , (58) 

L^I^^-ifM^ + l + r'a)]), (59) 

where \ represents a generic coefficient and the 0[s are the associated dimensionless coeffi- 
cients which govern the scale dependence of &[(/i): 

6^) = &rGuo)+AlrA (60) 

A 4 



Ecker et al.|38| have calculated all 0i in a theory without an explicit A resonance and 
obtained 

09 — 010 — 018 — 019 — 020 = 023 = 0, 

(61) 



0u = f + g\, 02i = g A + g%, 022 = -g\. 



(At first glance, the functions of |38| and are different but actually they are equiv- 
alent: Some terms which are proportional to the nucleon EOM are kept in ||39|| , but trans- 
formed away in [0q| . For convenience we adopt the functions of |38f| .) 

The A — 7r loop correction to 



«.</•) = *. " =^ + + Tab '« IT + (62) 



A = — + «/(— ) 2 " 1, (63) 



Now, collecting all other divergences, we get following results: 
G x (m) + 4(G 6 (/x) - G 4 (m)) = (?! + ^(G 6 - G 4 ) + A ?Y (9i- f + IS* ln 

(64) 



G 2 ( M ) + 4G 4 (m) = G 2 + 4G 4 + ^^f^(g A - \ 9l ){L + -L ln ^), (65) 

6 9 167r z /i 

where the Gi and G 2 are the coefficients associated with the 7AN vertex and the • over 
the G's indicates it is taken as lowest order value in the expansion — in other words, taken 
in the limit: m T — > 0, A — > 0, with — fixed. The remaining divergences are absorbed 
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by altering the The 7r° photoproduction amplitude can be renormalized only if the j3 
functions are modified in this way: 

2 (AttF^) 2 2 5 

Po + q M ^ 9nAN(Pc6 - Pa*) = ^91an(9a - ^9i) ■ (66) 

Pio = 0. (67) 

Here 

G l M = G l M+(3 Gi \n^. (68) 

/i 



The charged pion amplitude also requires an alteration of P21, P22 and P 
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(4vrF w ) 2 to a , 3 40 2 80 

#22 + 9nAN{PG6 ~ Pa) = ~9 A + ^fl^AiV ~ ^j^l&rAtf: < () 9 1 



fti = ^ - /3 22 , (70) 

a 9A i 3 50 2 2 2 

/?17 = y + #A - gjSl&rAJV + q9a9*an- ( 71 ) 



As Kambor pointed out [35], one cannot identify the coupling constants of theory including 
delta degrees of freedom with those in HBChPT, even if they multiply the same structures 
in effective Lagrangians. The reason is the process of integrating out the additional degrees 
of freedom leads to a (general infinite) renormalization of the bare coupling constants of the 
underlying theory So if we keep A finite, the bare coupling constants in our Lagrangians 
will differ with the usual one in HBChPT even in the chiral limit. From (83), (84) we have: 

Kv = K x_ m ^^ A - 2 ^ff gl AN {x\n.2x - Vx 2 - lm[z + Vx 2 - 1]), (72) 

K% = k v + 2 ^[ N gl AN (16n 2 L + In — ), (73) 

Gi = G\ + 2 7 5 2 m ;^ gi^A.(xln2x), (74) 

G% = G 1 - 2 §^ 9l w(16vr 2 L + In ^). (75) 
tin r w [i 

Here x=-^-. These constants ; uin and c\ are also infinitely renormalized: 

m N = m*- 4c*M 2 - MM + £m^1 R (A) 7 (76) 
1 * 32nF 2 \2tx 2 F 2 y MJ' y 1 



where 



R(x) = -A(x 2 - 1) 3/2 ln(z + Vx 2 - 1) + 4x 3 In 2x - 6x In x, (77) 

n 2 1 1 9 A 

m x = 7h + 90 A s { _ 16L + - In ^]), (78) 
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c? = ci + ^A(-16L + -li - In -]). (79) 



Another interesting example is mass of isobars |47]:(Here d\ is the coefficient associated with 
the counterterm : AATr(% + ).) 

m a = + A x - 4« - i^f - A), (80 ) 



864ttF 2 48tt 2 F 2 v M w 



2F 2 v 16vr 2 /i 8tt 2 " 
The renormalization of the axial coupling constant is similar but more complicated j48 



4m 2 h r A 

2a = <M + 27?r2F2 A {16tt L + In — + -}, 



where 



S(x) = ai + a 2 x 1 + a 3 x 2 \n2x + a^xy/x 2 — 1 ln[x + \/ x 2 — 1] 

+a 5 I ^^ln[x+v / ^I], 
-4 3 136 2 200 2 

a l — "3"^ ^~9Ag n AN + T^fi'l^TrAiV' 

_ 647T 2 
a 2 - -^-gA9irAN> 

304 2 400 2 

a 3 ^j-9Ag n AN + -^9l9nAN, 

Ifi 2 400 2 

128 2 

a 5 ^j-9A9 n AN- 



81) 
82) 

83) 
84) 

85) 

86) 
87) 
88) 
89) 
90) 



In general, the LEC's in the present expansion are different from HBChPT without an 
explicit A degree of freedom. The reason for this is obvious: processes including the explicit 
A contribute and serve to renormalize the parameters. However, in the chiral limit of m n = 0, 
k%, G\, m x , Ci and g\ can be identified with the analogous LECs in HBChPT without A's. 
In contrast, parameters of k v , Gi, c/a, ttin an d ci are defined in the limit: A — > 0, m n — > 0, 
-A fixed. 

There is an interesting formal limit to consider apart from the chiral limit; namely the 
large N c limit. Recall in this limit one naturally has A — > with m n finite. Moreover, in the 
large N c limit, qa = §#1 = ^p<?7rAAT- Inserting this into our expressions for the /?'s we find: 

At = y- (91) 



14 



Actually such a result is required by the consistency of N c counting rules. Note that qa 
is 0(N C ) quantity, and F n ~ 0{^/W c ) and m n ~ 0(1); therefore all terms with g\ must 
vanish in the large N c limit. Such a cancelation is possible only when the A-isobar degrees 
of freedom are included |51t],|52|. 



6 Comparison with Experiment 

The values of unknown parameters such as G\, G^--- are expected to be extracted from the 
experimental data, and to be used to make predictions in other processes. However, as 
mentioned in the Introduction, there is an obstacle: The amplitudes at the A pole diverge, 
and amplitudes in the vicinity of the pole clearly cannot be taken seriously. 

At first glance, this appears to be a fatal problem for this approach since our interest is 
precisely in the A region. One natural way to cure this is to use a dressed delta propagator, 
i.e, put the self-energy part T,(E) in the delta propagator. With a self-energy included, 
amplitudes become smooth since the imaginary part in T,(E) shifts the pole from the real 
axis. However this approach breaks the power counting scheme in its purest form because it 
requires that part of the interaction (those associated with A decay into pion plus nucleon) be 
iterated to all orders while other parts are not. The overall expansion becomes questionable, 
and as we pointed out in the previous section, the power counting scheme is our only way 
to make predictions consistently. 

At a technical level, the problem is that away from the A pole, the A propagator is (9(e _1 ), 
while in the immediate vicinity the pole is 0(e~ 2 ). As one moves up the resonance, the 
behavior changes, making a systematic treatment problematic. One purely phenomenological 
alternative is to simply insert the empirical decay width F in the invariant amplitude as 



Adelseck et al. have done[56j. Such a scheme is not systematic, however; and moreover, 



unitarity is violated 45[ . 



Therefore we suggest that instead of directly comparing the theoretically calculated am- 
plitudes with the experimentally extracted ones at all energies, we only compare weighted 
integrals of the amplitudes. For example, 

M$ = / M l± (E)W n (E)dE. (92) 

Here E is the energy of the photon in the cm frame, W n (E) is a smooth weight function, 
and n is an integer index specifying the particular choice of weight function. 

To justify such an approach, the weight functions need to be chosen with some care. First 
of all, the hierarchy introduced by power counting must be maintained after integration. In 
other words, the amplitude with factor E _^ +ie should not be enhanced beyond what is 
permitted in our power counting. To satisfy this requirement we must integrate through the 
entire A region. Actually there's another reason to do so: Recall that in the vicinity of the 
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pole, the theoretical resonant amplitude has a 5 function in the strength through 



lim \ = P(— L— ) - m8(E - A), 

but the actual resonant amplitudes will be more like a Breit-Wigner form and one expects 
the imaginary part to spread the width T. Therefore, the weight functions are required to 
cover the whole A region to keep the full information of the imaginary part in an experi- 
mental resonant amplitude. On the other hand, E max could not be put too high. The C(e 4 ) 
contribution of amplitudes becomes more important when energy increases and our calcula- 
tion loses its predictive power in the higher energy region. To satisfy both requirements, the 
best place for E max is at the upper end of the A region. Since we are interested in power 
counting it is sensible to look at a Taylor series for the W n (E): 

W n {E) = 0?) + ,£> A + a ( 2 )(A )2 + a (3) ( A)3 + a (4)( A)4 + ... (93) 

/ / Lp 1 1 Lp 1 1 bp 1 1 bp 

It is clear that the effect of higher order coefficients such as do not contain reliable 
information since they can not be distinguished from the effect of 0(e 4 ) amplitudes which 
are absent from our calculation. Accordingly we can choose 

E E E 

),(— ) 2 ,(— ) 3 , 
m p m p m p 

as our "basis functions". Any W n (E) is equivalent to their linear combination once the 
higher order terms are thrown away. Thus we choose W n (E) as 

W n (E) = (— )». 
m p 

The preceding analysis also suggests another advantage of studying these integrated quan- 
tities apart from the problem of the A pole. Note there is a finite amount of information 
which can be extracted from a systematic expansion at a finite order. Direct predictions 
for cross sections as functions of energy formally have an infinite information content since 
it takes an infinite number of parameters to describe an arbitrary function. Clearly, much 
of the information content contained in a predicted functional dependence has considerable 
correlations. It is useful, therefore, to construct a number of discrete observables, such as 
our integrals which characterize the energy dependences. The scheme proposed here corre- 
sponds to picking the maximum number of independent predictions which we can make at 
this order. 

We integrate from the threshold through the entire A region: 

_ /„\ I rEmax Ej 

M[ n J = — / M l± (E)(— ) n dw. (94) 
The additional factor l/m p is only to ensure that M. has the same dimension as M.. The 
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values of the unknown parameters can be fit from the amplitudes as follows: 
RePf p = eg A (f + eg A k p (* B + eg A K p ^ R + eg A (l + «p)ciCf 
+ e 9\CP + e 9nANGiC^B + e 97vANGiC^R + eg 7T ANG 2 Cf + eg^^G^ 

(95) 

+eg nAN G 1 C l H + egl AN g A (? + egl AN9l tf + efe 9 Cf , 
i = 1,2,3. 

\lmPf p = ecfi&P + eg nAN G^ B + eg^ N G^ R + eg vAN G 2 Z? + eg^ AN G£f 

+eg^ AN G^ + egl AN g A £? + eg 2 wANgi tf, (96) 
i = 1,2,3. 

Here G 2 = G 2 + 4G 4 , G 6 = G 6 - G 4 , c x = m N c 1 , k p = \(k v + K 8 ),b g — b 9 — b w - 
^Irr? 9k&nGi(1 + Ax + Az + Ylxz). k p means the parameter taken in the limit: A — > , 
m ir - ► 0, fixed, and k p ; the 7r loop correction is included in this. 

The first four terms in (146) are from tree graphs without delta; the sixth to eleventh 
terms are due to tree graphs with delta. Note that such tree graphs also contribute to the 
imaginary parts of amplitudes due to the delta function in E _ A+ie - The fifth term is from 
loop graphs without delta; the twelfth and thirteenth terms are A — ir loop contributions; 
the last term, which only appears in P 3 , is due to the counterterms in C^ N . Note that 
the quantities (as they must be), like are //-dependent, however final amplitudes are 

independent of \i because the k v , G±, G 2 and G 6 are also /i-dependent, and compensate the 
ones from the loop. Note that the contribution of G e always couples with the one of g wA N, 
it's impossible to separate them by fitting the pion photoproduction data. But only the 
latter appears in the n — N scattering amplitudes, therefore their values can be determined 
by both experimental data. 

It is obvious that power counting is preserved in this scheme with this set of weight 
functions. The 0(e 2 ) amplitudes represented by (f, (? R or (^ R , (i=l,2) are significantly 
larger than the other terms, and or (^n an d ([ ', C G are n °t particularly enhanced which 
shows that the wild behavior of E _ A+ie is tamed by our weight functions. 

We set A=294 Mev, F ff =92.4 Mev, M^=938 Mev, E max =3A0 Mev, and ^=500 Mev. 
The following results are given (in the unit of 10~ A /m n ): 
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n 






r B 

Sl;B 


Ci 


Ci> 




Sl;B 


cr 


Ci 


C? 


C K 


S1;A 


C L 

S1;tt 


S1;A 





23.89 


26.93 


-1.04 


2.40. 


9.23 


-12.07 


-3.92 


0.74 


-2.96 


-2.96 


-1.07 


15.62 


6.75 





1 


6.60 


7.52 


-0.36 


0.67 


2.66 


-4.65 


1.03 


0.29 


-1.18 


-1.18 


-0.72 


4.31 


3.04 





2 


1.90 


2.19 


-0.12 


0.20 


0.80 


-1.67 


-0.28 


0.11 


-0.44 


-0.44 


-0.31 


1.18 


1.17 





3 


0.57 


0.66 


-0.04 


0.06 


0.24 


-0.58 


-0.08 


0.04 


-0.16 


-0.16 


-0.12 


0.32 


0.42 





n 


Zi 


VL;i? 


S1;B 


£l 


ii 




CK 

Sl;ir 


cK 

S1;A 


f L 

S1;tt 


£l;A 







-3.47 


6.86 


-0.95 


-0.47 


1.90 


1.90 


2.51 


2.53 





-6.53 


1 


-0.97 


1.90 


-0.26 


-0.13 


0.52 


0.52 


0.98 


0.70 





-1.80 


2 


-0.28 


0.52 


-0.07 


-0.03 


0.14 


0.14 


0.35 


0.19 





-0.50 


3 


-0.08 


0.14 


-0.02 


-0.01 


0.04 


0.04 


0.12 


0.05 





-0.14 



n 


<>2 


^2;R 


^2:B 


^2 


C D 

S2 


C E 

S2;R 


C E 
S2;B 


S2 


<>2 


C H 

S2 


t >2; 7 r 


<*2;A 


t >2; 7 r 


C^A 





-41.53 


-24.82 


1.66 


-2.40 


-12.10 


12.01 


5.48 


1.49 


2.96 


2.96 


0.54 


5.82 


-5.36 





1 


-11.41 


-6.91 


0.54 


-0.67 


-3.46 


4.65 


1.93 


0.59 


1.18 


1.18 


-0.31 


1.61 


-1.94 





2 


-3.27 


-2.00 


0.18 


-0.20 


-1.02 


1.67 


0.66 


0.21 


0.44 


0.44 


-0.21 


0.44 


-0.67 





3 


-0.97 


-0.60 


0.06 


-0.06 


-0.31 


0.58 


0.22 


0.07 


0.16 


0.16 


-0.09 


0.12 


-0.22 





n 




S2;il 


tE 

S2:B 


$ 






cK 

S2;tt 


cK 

S2;A 


t L 

S2;tt 


t L 

S2;A 







0.80 


-6.86 


-1.89 


-0.95 


-1.90 


-1.90 


0.07 


2.17 





2.56 


1 


0.23 


-1.90 


-0.52 


-0.26 


-0.52 


-0.52 


0.12 


0.70 





0.71 


2 


0.07 


-0.52 


-0.14 


-0.07 


-0.14 


-0.14 


0.06 


0.17 





0.19 


3 


0.02 


-0.14 


-0.04 


-0.02 


-0.04 


-0.04 


0.02 


0.05 





0.05 



n 


ci 


Cs% 




C 3 G 








c 3 G 




ct 


^>3;A 




S>3;A 







-0.36 


3.68 


-7.50 


-9.32 


12.41 


-1.88 


1.16 


9.28 


9.28 


-0.67 


9.80 


4.69 





19.26 


1 


-0.16 


1.06 


-2.23 


-2.50 


5.94 


-0.41 


0.41 


3.28 


3.28 


-0.19 


2.70 


2.90 





5.59 


2 


-0.06 


0.32 


-0.68 


-0.70 


2.35 


-0.09 


0.15 


1.20 


1.20 


-0.05 


0.75 


1.24 





1.68 


3 


-0.02 


0.10 


-0.21 


-0.20 


0.86 


-0.02 


0.04 


0.16 


0.16 


0.02 


0.22 


0.47 





0.52 


n 


f E 








£3 






tK 
S3;A 


?3;tt 




;A 







-13.72 


-0.95 


-0.47 


-3.80 


-3.80 


0.90 


0.36 





-9.09 


1 


-3.79 


-0.26 


-0.13 


-1.04 


-1.04 


0.42 


0.10 





-2.50 


2 


-1.04 


-0.07 


-0.04 


-0.29 


-0.29 


0.17 


0.03 





-0.69 


3 


-0.29 


-0.02 


-0.01 


-0.08 


-0.08 


0.06 


0.10 





-0.19 



Table 2: £{? A , £f, t;f and represent the contributions due to the delta function of E _ A+ie in 
the loop diagrams like (2-A-l). The other £j's are from n loops. 

Unfortunately, at present there is not enough data to test this theory. This scheme 
only predicts integrals of amplitudes through the Delta region. In order to compare with 
experiment it is essential that the amplitudes be separated out so that these integrals can 
be estimated. To do this, more precise and complete measures of spin observables are 
required. The cross section plus three single-spin observables determine the magnitudes 
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of the amplitudes, but double-spin observables determine their relative phase. By carefully 
selecting four of the double-spin observables, one can extract all of the requisite phase without 



discrete ambiguities [50, ],[55|,[53|. Once the experimental data of amplitudes are available, 



the values of these unknown parameters can be extracted and we can test our predictions. 

A fundamental issue is the predictive power of our method of comparison. The scheme is 
designed to separate out the maximum number of independent quantities extractable from 
experiment at a given in the small scale expansion. To be predictive, there must be more 
observables than free parameters. In neutral pion photoproduction, the P-wave amplitudes 
are determined by eleven unknown parameters. With these eleven parameters we have 12 
integrated observables to fit. Fortunately from the S-wave amplitudes: 

Ef + P (uj) = E^ n ' N + E^- nLoop + E*° TnA . (97) 

2 4 
^Born,N e 9A , eg A r „ 2 , a 2 771 tt n a~ \ , .• (a, .2 2' 



£o+' JV = ~ 9A F i— + M + J3 2 F {6^ + 4 < - " 8 £ i) + *p(^ - <))■ 

247im N F n uj 48iim A N F n u z 

(98) 

E 0+ P = Q^2 F 3 K arCSil1 — - U \l m l ■ (") 

^BornA _ e&rAivGi cj 5^ 2 - 2m^ ^(u; 2 - ml) 
Eo+ -4vrm 2 v F^ + A [ 27 + 9(u + A) J ' 

here uj is the pion energy in the cm frame. Four parameters: gA, c\ and g n ANGi can be 



extracted without integrals because that our calculated E§+ doesn't diverge when u — > A, 
therefore only 7 unknown parameters need to be fit from P-wave multipoles. Toward the 
charge pion photoproduction, the P-wave has only one new parameters: 622 + ^23, and again 
theoretical results of E^^ never diverge in A region therefore four parameters: g r A , 617 — 
621 and k n can be determined by matching them with experimental data. 



Furthermore if we generalize to the case of electroproduction, there are only two addi- 
tional parameters which need to be fit, and the number of observables increases since there 
are additional C2 amplitudes. Therefore up to 0(e 3 ) the predictive power of our approach is 
clear. If one works at higher order, the number of basis functions increases, therefore a larger 
number of independent observables becomes available. However many new parameters will 
be involved. Whether one ultimately has increased the predictive power has not yet been 
settled. 

In conclusion, the "small scale expansion" provides us a systematic way to calculate the 
processes of the nAN system, and with the power counting scheme we can fit the unknown 
parameters to make predictions. Our method is designed to isolate the maximum number 
of independent predictions which can be made at a given order. 
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Appendix A: Feynman Rules 



Here we collect the Feynman rules which are needed to calculate tree and loop diagrams. 
The following notations are used: 

1: momentum of a pion or nucleon or delta propagator; 

k: momentum of photon; 

q: momentum of external pion; 

e: Photon polarization vector; 

p: momentum of a nucleon in heavy mass formulation; 
r: momentum of a delta in heavy mass formulation. 

Isospin indices of pion are a, b, c, d..., the isospin indices of A are i, j, k..., and the spin 
indices of A are /i, z/, a. is the nucleon four-velocity and is a covariant spin- vector. We 
also give the orientation of momenta at the vertices, i.e, which are "in"-going or "out" -going. 
Q=diag(l,0). Here we only present the ones related to A since others can be found in 

Vertices from £^ : 

(1) A propagator: 



1-3 

fiu S,/=3/2 



(2) One pion (q out): 

(3) One pion, one photon: 



v ■ I - A + ie 



-^€L V 2T%t V2 {S. q )g : 

■F TT 



fJLV 



- %e 9l cSbtik bpkj (q,\ f . 
p 6 S/=3/2 r ^/=3/2l D € )9^ 



TT 



Vertices from 

(1) A propagator: 



(v-l-A + ie) 2 2m N 



Vertices from C^ N : 
(1) One pion (q,out): 



(2) One photon (k,in): 

(3) One pion, one photon: 



ted 



m N 



d^KS-k^-iS-e)^ 



ieg^AN c3i 

Fir 
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Vertices from C^ N : 
(1) One pion (p in; r, q out): 



9wAN i U d QirAN i uct 

(2) One photon (k,p in; r out): 
" ( u . A; )^[-/ ll (5. e )^ + ^(5-r) e ^-^ 2 (5-r) e ^V«^ + ^(5-r-p)-G ieM (5-r+p)] 



2///-;. ' L '~ iV ~ ~ r *' 2 ' ' '~ '■ " ' 2 



it 



+^(e-v)6 i3 [h 1 (S-k)r„-f(S-r)kv+ih 2 (S-r)e^ 



+— (v ■ k)S iS [G 6 (S ■ k)e, + (G 4 - G 6 )(S • e)fc M ] 
+— (e • v)5 i3 [-G 6 (S ■ e)k^ - (G 4 - G 6 )(S • fc)e M ] 



foi = -2d - 2xoGi + y C2; ^2 = yoG 2 - 2x Gi. 



(3) One pion, one photon (p,k in; r,q out): 



r9nAN -e c3 \q,(e-v) + e,(vq)). 
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